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Abstract 

The low lying spectrum of QCD in the (5-regime is calculated here in 
chiral perturbation theory up to NNL order. The spectrum has a simple 
form in terms of the pion decay constant F and a combination of the low 
energy constants Ai and A2. Since measuring low lying stable masses is 
among the easiest numerical tasks, the results should help fixing these 
parameters to good precision. 



1 Introduction and summary 



The low lying spectrum of QCD, in a special environment, is that of a simple 
quantum mechanical rotator [1]. Although amazing as it is, this spectrum will 
never be measured in QCD experiments. This physics, however can be studied 
numerically which will deliver precise predictions on some of the low energy 
constants in chiral perturbation theory [21 131 H] . One of the reasons is that 
measuring low lying stable masses is among the easiest numerical tasks. Further, 
the QCD rotator lives in a box of size Lg x Lg x Lg x {Lt oo), which creates 
an infrared-safe environment^. This allows to study the chiral limit first and 
switching on the symmetry breaking terms later. In two-flavor QCD, in the 
leading order (L) of chiral perturbation theory the SU{2) x SU{2) ^ 0(4) 
rotator has an inertia 9 proportional to the size of the spatial box: 9 = F^L^. 
Here F is the pion decay constant in the chiral limit. In the next-to-leading order 
(NL) the inertia is corrected 9 = F^L^g{l+ - l/F^L^) ff], where l/F^L^ ig 
the small expansion parameter. The logarithms and the additional low energy 
constants li and I2 enter only in the NNL order, which is the highest order we 
consider in this work. 



In the environment discussed above (J-regime [T]) it is natural to divide the 
degrees of freedom into fast and slow modes. The fast modes can be treated 
in perturbation theory, while the slow modes build the slowly moving rotator, 
whose energy excitations are much smaller than those of the standard Goldstone 
boson excitations which carry finite momenta. 

Let us summarize the final results before going over the details. We quote 
the result for N = 4, which corresponds to 2-flavor QCD. Up to NNL order 
the rotator spectrum in the chiral limit has the form □ 

Ei^^l{l + 2), 1 = 0,1,2,... , (1) 

where N refers to the underlying group 0{N). The corrections from the pcrtur- 
bative expansion enter in the inertia 9 which, up to NNL order in the chiral 
limit, reads 

9 = F'^Ll I 1 - — J-^G* 



1 



(F2L2)2 



0.088431628 (2) 
1 /I 



^ dOdOG*^ (^^ ln(AiL,)2 + ln(A2L.)2 

^It is known since a long time in condensed matter physics also that in a finite box at zero 
temperature the lowest excitations in the spontaneously broken phase are related to the slow 
precession of the order parameter described by a rotator ref. |5l|6]. 

^Terms with other Casimir forms are expected to enter beyond NNL order. 
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Here Ai , A2 are the standard scales related to the bare low energy constants 
li, I2 in the chiral Lagrangian. The conventions for Ai, A2 are given in section 7. 

We are using dimensional regularization (DR) in this work. The constants 
G* and dOdOG* are related to the constrained Green's function D*{x) and its 
second time derivative dodoD*{x), which enter the perturbation theory: 

D*{0) = -^G* , G* = -0.2257849591 . (3) 



dodoD*{0) = -^dOdOG* , dOdOG* = -0.8375369106 . (4) 

In the Green's function D*{x) the non-pcrturbative slow modes (rotator 
modes) are missing (notation D*) and the UV singularity is also subtracted 
(notation D). The precision of the numerical numbers are estimated to be 
10~^, or bettejfl. Definitions and properties of the Green's functions we use are 
summarized in section 5. 

Although the final result in eq. ^ is very simple, the underlying chiral 
perturbation theory is not. For this reason it is a good news that an independent 
calculation is under way using a completely different technique [5]. Including 
the symmetry breaking contributions is also in progress [S]. 



2 The chiral action and the rotator in leading 
order 



The low energy limit of QCD with two massless quarks m„ = = is de- 
scribed by an effective non-linear 0{N = 4) tr- model. The Lagrangian, up to 
1-loop level, has the form L^s ~ L'^g + , wher^ 

- -h (6)^s9^s)(a,sa,s) - h (a^sa,s)(9^sa,s) . 

Here F, li , I2 are the bare low energy constants and the A^-component field has 
unit length 

Saix), a = 0, 1, . . . , (A^ - 1), S2(x) = 1 . (6) 

•^Of course, no such precision is needed in this perturbation theory. It might help though 
when comparing with related works in the future. 

*We assume that rotation symmetry is respected by the regularization. In dimensional 
regularization this is the case. 
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The action in eq. ([5]) might also be interpreted as the effective low energy 
prescription of a ferromagnct. We consider rh Lg X Lg X Lg x Lt box with 
spatial volume Vs — L^, while the (Euclidean) time extension is taken very 
large, Lf — oo. In the — > oo limit the system has a net magnetization and 
massless Goldstone bosons. 

We consider a cylinder geometry ((5-regime), where Lg is finite, but suffi- 
ciently large so that the (would be) Goldstone bosons dominate the finite size 
effects. Due to the microscopic magnetic moments, the Lg x Lg x Lg spatial 
box has a net magnetization on each time slices. Since Lg is finite, this net 
magnetization is moving around as the function of the time t. 

In leading order (L), which is a good approximation if the dimensionless 
expansion parameter 1/{F^L'^) is small, the microscopic magnets on a time 
slice are parallel and the rotator action from eq. © goes over to 

A^t^ ^ J dte{t)e{t). (7) 

Here e{t) is the direction of the total magnetization in the internal 0{N) space 
at the time t: 

e(<) = [e(t)o,e(t)i,...,e(i)jv-i], e{tf = 1 . (8) 

Equation ([7]) describes a quantum mechanical 0{N) rotator, the QCD rotator 
in leading order [T], with a discrete energy spectrum above the ground state: 

U=^ , e = F^Vg, Ei = ^lil + N-2), / = 01 ,2 ... . (9) 

In eq (ini) H is the Hamiltonian, L is the 0{N) angular momentum and Q is 
the inertia of the rotator. Our aim is to determine the corrections up to NNL 
order, where the low energy constants ^i, I2 first enter. 



3 Separating the slow and fast modes 
3.1 New integration variables 

If Lt ~ , it is sufficient to take special care of the freely rotating total magne- 
tization [ini HH HH [13] . For Lt/Lg large, however the magnetization on distant 
time slices in the cylinder might differ significantly. We have to treat these 
special slow modes non-perturbatively. These modes are the k = (/cojk = 0) 
modes in Fourier space. In this subsection we deal mainly with the measure, 
while the action is treated in Sections 3.2 and 3.4. 
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The steps followed here are similar to that used in [7], but it is simpler and 
exact in every order. Unlike in [7], where lattice regularization was used, we 
apply here dimensional regularization. 

Insert into the path integral the identity 

1 = n / dmit) n '5K(t) - - / S-{t,^)] , (10) 

t a=0 " 

where x is the spatial coordinate, x = (i,x), and 

m(f) = m{t) e{t) , = 1 , dm(t) = ■m{t)^^^dm{t) de{t) . (11) 

The 0{N) vector of unit length e(<) is the direction of the 'magnetization' 
J^S{t,x.) on the time slice t. The local 'magnets' S(t,x) fluctuate around the 
slow mode e{t). 

Introduce the 0{N) rotation i7(t) as 

e(t) = r!(i)n, n=(l,0,...,0) . (12) 

Having an 0{N) matrix D, rather than a vector will be convenient in the ma- 
nipulations below. We consider as a function of e{t). Note, however, that the 
path integral 



Z=YlJ dS{x)SiS'^{x)-l)Y[ J dm{t)m{t)^-^ j de{t) 
5^[m{t)Vt{t)n~ — f S(t,x)]e-^'=«(^) , 



(13) 



depends only on the first column of the matrix fl. 

Introduce new integration variables R in the path integral: 

s(t,x) = n{t)j:{tfR{t,x.) , (14) 

where E(t) e 0{N) and R has unit length. The matrix I]{t) is taken to have a 
special structure: S(t)y — S(i)ij, ^{t)oo = 1, ^{t)oi = 0, S(i)io = 0, where T,{t) 
is an 0{N — 1) matrix. This 0{N — 1) matrix E(i) will be chosen conveniently 
later. 

The partition function reads 

Z=Y[ f dR.{x)S{R\x)-l)Y[ I dm{t)m{tf~^ f de{t) 

* (15) 
5^[m(i)r!(t)n - n{t)j:{tf^ f R(t,x)] e-^e«(f2S^R) 
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Using 5^ {riz) = (z) and En = n gives 

Z^YlJ dRix)S{K^{x)-l)Y[ J dm{t)m{tf-'^ J de{t) 
6^[mit)n ~ - [ R(i,x)] e-^o«(os-R) 



(16) 



eq. shows that on each time shoe the vector R fluctuates around the 

(1,0,..., 0) internal direction: 

R(i, x) = (1, 0, . . . , 0) + small fluctuations . (17) 

The small fluctuations can be treated in perturbation theory: 

R(t,x) = [(i-n2(i,x))5,n(f,x)], (18) 

where 

n(t,x) = [n(t,x)i,...,n(t,x)^v-i] (i9) 

are the small fluctuations. From eqs. ([1^)) . follows 

/ n(t,x), =0, t = l,...,N-l. (20) 

The field n(a;) is small. This fact allows to integrate out the Il-fields ('fast 
modes') in a systematic perturbation theory. We should keep in mind that the 
slow modes e{t) were separated. The k = (A;o,k = 0) modes in Fourier-space 
are slow and arc not part of the fast Il-fields. 



Two of the intergrals in eq. ()16|) are zero in dimensional regularization. The 
first one is the well known measure which is created as we replace R by the II 
integration variable. The second integral, when put in the exponent, has the 
form (N - 1) ln[l/Vs /^((l - U'^it,x))'^^^)], which is zero also. The relevant 
part of the partition function Z will be given in Section 4. 



3.2 The slow and fast modes in the action ^^^(fiE^R) 

Here and in Sec. 3.4 we consider the action which is a classical object. In 
the following steps we found it useful to introduce an infinitesimal parameter e 
which disappears at the end. 

In Sec. 3.1 we introduced new variables S(i,x) — i7I]^R(i,x). Let us write 

(2) 

the time derivatives in the leading action A^^^ in the form 

9tS(i,x)9tS(i,x) = hl-S{t + e, x)S(t, x)],^n- (21) 
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Eqs. (mi), dm) imply 
9^S(i, x)9^S(t, x) = a^R(i, x)9^R(i, x) - ^Q{t)R{t + e, x) R(t, x)e^o , (22) 

where the N x N matrix Q{t) is expressed in terms of fl and S 

Q(<) = y(< + e)S(t + e)^ - 1. (23) 

Here we introduced the notation 

V{t + e) = J:{t)fl{t)'^n{t + e) . (24) 

The 0{N) matrix I](i) is defined in terms of the 0{N — 1) matrix S(t) (see 
after eq. (fll)) ). We fix this matrix now: 

It is easy to show that S G 0(A^ — 1), indeed. Further, Q{t) can be expressed 
now in terms of V{t + e): 



Q{t)oa = V{t + e)oo - 1 , Q{t)o, = -V{t + e)^o , Q{t)^o = V{t + e),o , 



V{t + eUV{t + e),, (26) 



1 + V{t + e] 



00 



Assume now that V is known at some tp, while fl(t) is known for any 
Having V{to), we can fix S(<o) and S(io) using eq. P5|) . Actually, we have the 
following chain 

^(io) S(to) -^^(io + e) -^S(to + e) Q(to) ■ (27) 

eg. |25j eq. eq. l25j eq. J26j 

At the and of this chain we have Q{to). In addition, we have V(tQ + e), so we 
start the chain again creating Q(to + and so on. 

We shall see that for the rotator spectrum up to NNL order the following 
three combinations of the matrix Q are needed only: Q{t)QQ , Q{t)ioQ{t)io and 
Q(t)i,i (the repeated index i is summed). These combinations can be expressed 
in terms of the slow modes e{t) in the e ^ limit: 

\Q{t)m - -\e{t)e{t) , \Q{t),,, - -;^e(t)e(t) , 

'I 2 62 2 (28) 

-Q(t),oQ(t).o = e(t)e(t) . 
Let us demonstrate eq. 



'It will soon turn out that only the first column (i.e. the vector e) is needed. 
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00 



Using eqs. ([26|) . (|23l) and So, a = (^o^a we obtain 

^QWoQ = + e) - l)oo = + - l)oo • (29) 



e2" 



Expanding in e gives 



^[f7(t)^f}(<)]oo = ^^mfm]oo = -le{tm . (so) 

Here we used 

2n{t)'^n{t) + n{tfn{t) + n{t)'^n{t) = o . (3i) 

Turn to the next case: 
Consider first 

-Qit)ao^-[^{tMtfnit + e)-^{t)U^[^it)n'^{t)]abmb, (32) 

e e 

where we used that ^{t)ao = ^ao- It follows then 

-Q{tU-Q{t)^o = -Q{t)ao-Q{t)ao = e(t)e(<) . (33) 
e e e e 



Here we used that ^Q(i)oo is 0(e). 
eq. p6)) implies 



Here F(t + e)iQ — Q{t)iQ and the denominator is 2 + O(e^), we obtain from 
eq. ([55)1 the result eq. 



As discussed before, the fast modes are carried by the R = [(1 — n-^)5,n] 
field, where the k = {kg, k = 0) modes are missing. The first term on the r.h.s. 
of eq. ((55]) gives the fast H — H interactions, while the second term describes the 
fast-low interactions. We get for the leading action A^^g in eq. ([5]): 



El 

2 2 (35) 

^g(t)oo[i - u'ix)] - -^Q{t\,n{x\ \i{x)j 

2 
-( 

e 



-[d^-R{x)d;R{x) 



'-Q{tu[dt{i - u\x)Y^ n{x), - dtn{x\ (i - u\x))^ } 
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The terms above give the action A^^ff up to NNL order. We shall expand in the 
Il-fields in perturbation theory. The last term in eq. (|35|). which is odd in 11, 
enters only in the NNL order of this expansion. 



3.3 Counting rules 

The small expansion parameter in the ^-regime is 1/F'^LI = 0{5^). The expan- 
sion of the rotator action has the form 

^VMmm+ - ^+ - jpj^y, ■■■) (36) 

in the leading L, next-to-leading NL, NNL,. . . order. 

The finite part of the pairings of the fast modes are < HH > ^ 1/{F'^L'^) 
and < d^Ud^U > ^ 

In the expansion there are terms also with quadratic and higher powers of 
the slow mode e. The following consideration shows that an additional ee 
term in the bracket of eq. ([55)) is 0{S^), i.e. beyond our NNL calculation. 

The argument is as follows. The leading Lagrangian is ^ F^L^ee which gives 
for the conjugate momentum L ~ F^L^e. It follows then ee ~ h^{F^Ll)~'^. 
Since 0{1), we obtain 

e{t)e{t)Ll^i^r^O{S*), (37) 

giving ee ~ 0{6^) as stated above. In the NNNL order, which is beyond our 
calculation, such corrections are expected to enter. Consider an example. As 
we shall see later, in the Boltzmann factor enters (among others) 



^ ^e(i)e(i)n(x)n(a;)] = 1 + ^ ^e(t)e(t)n(a;)n(x) 



exp 

J X 

1 r F^ r F^ 



2 /., 2 



■e(t)e(<)n(x)n(x) f —e{t')e{t')U{x')U{x') 
Jx' 2 



(38) 



We pair out the fast modes on the r.h.s. of eq. (|38l) . Pairing n(a;) with n(a;) and 
n(a:;') with n(a;') in the second line is just needed for the exponentialization. 
Pairing n(x) with n(x') gives, however 



[mmmmn / / D*{x~x'f = 

' (39) 
jj^^^ F'Ll [emt)] [e(t - t")e(t - t")] £ D*{x"rj^ . 
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In D*{x") the k = is missing, the smallest k is ~ l/Lg. Therefore, D*{x") has 
an exponential cut on the level t — t" ~ Lg. The time distance is small for the 
rotator, therefore e{t — t") ^ e{t). The integral over D*, after renormalization, 
is 0(1) giving 



3.4 The action A^;J(^]S^R 



In the term with li in eq. we can use the result in eq. 



d^Six)d^S{x) d,S{x)d,S{x) - a^R(a:)9^R(a:) d,Rix)d,R{x) 
2 

2 df,R{x)df,Rix) ^g(t)R(t + e, x) R(t, x) + 
^Q{t)R{t + e, x)R{t, x) ^Q{t)R{t + e, x)R(t, x),^o ■ 



(41) 



The leading term of the first part on the r.h.s. in eq. (PT|) is 9^119^11 9^119^11. 
This term can influence the low mode spectrum only if it is multiplied with 
both low and fast modes like eellll. This combination is is far beyond our 
calculation. Up to NNL order we can write for the li part 

d^S{x)d^S{x) d,Six)d,S{x) = ~^Q{t)oo d^U{x)d^U{x) 

4' (42) 

+ ^Qit)ro Qit),odon{x),doii{x)j . 

The term with I2 in eq. ([5]) requests somewhat more, but trivial work. It has 
the form 



d^Six)d,Six) 9^S(x)9,S(2;) 



1 ^ 



(43) 



a=l 



The part of the action which is relevant up to NNL order has the form 
Ai^m^R)^ 

^Q{t)oo 9on(x) 9on(a;) - ^Q(t),o Q(Ojo5on(x), 9on(a;), 



+ I2 



1 ^ 

QitUQit)jo{dQTl{xhdon{x)j + -J2daTi{x),da,nix), 



Qit)oo doU{x) doUix)]} 



(44) 
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4 The path integral up to NNL order 



The relevant part of the action is the sum of the terms on the r.h.s. of the equa- 
tions and (j33|. This action goes in the path integral as exp(— A^^^ —A^^), 
where the action depends on the slow e(t) and the fast n(x) degrees of free- 
dom. We consider the partition function eq. ()16p in dimensional regularization 
(DK,MS). Among other conveniences, this simphfies the measure. The Boltz- 
mann factor exp(—Al^g — A^^) is expanded in the 11 fields: only the leading 
parts remain in the exponent, the rest is a Taylor expansion in powers of 11. 

The path integral has the following form up to NNL ordei0: 

Z = 

N-l 



(1) 

(2) 
(3) 

(4) 
(5) 



exp [- ^ ^VsB{t)e{t)\ exp [- ^ ^a^n(x) d^Jl( 
[l + j^^e{t)B{t)U\x) + ...] 



x) 



F^2 



,1 



F'-QitU{-n'ix)doU{x),-U{x)don{xMx), 



F^-^Q{yo)jo{ln\y)doIl{y), - Uiy)doU{y)U{y),)] + . . . } 



(45) 



{Uix)d^U{x)){U{x)d^Uix)) 



{l + 4 ?i ^ hit)e{t)d^Ilix) 9^n(x) 



} 



(6) { 



+ ^QitUQit)jodoUix),doU{x)j + 
le{t)e{t)doU{x) daUix) 



e 

l+4?2 



} 



1 ^ 

{doii{x), doii{x), + - ^ d^ii{x),dji{x)j^ + . . . } 



The integration variables in this path integral are the slow modes e{t) and 
the fast n(a;) modes with the constraint 1/14 n(i, x)^ = 0, i.e. in Fourier 

^To make sure the notation: the brackets . . . > are multiplied with each other. 
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space the (fco,k = 0)-modes are missing. The corresponding pairings are 

< n(a;),n(0), >= ^-^D*{x) , < dodoll{x),n{0), >= ^ dodoD*{x) . (46) 



of 



In eq. (|45j) the six hnes numbered as (1) • • ■ (6) correspond to the expansion 
the action in the exponent. We gave exphcitly the terms only which are 
needed up to NNL. We integrate out in eq. the H-fields and obtain a path 
integral in quantum mechanics for a rotator. The first two contributions below 
are the other five are NNL: 

(1) 

[l + J^lF'VMt)e{t)^D*iO)} 
(2) 

(1) — (4) crossing 

{l - '^F^VMt)e{t)^D*{0)D*{0)} 

(2) — (4) crossing 

{l + J^lF^VMt)e{t)^D*{0)D*{Q)} 
(3) 

{l- J^lF^VMt)e{t)^ LdodoD*{z)D*{z)D*{z)} 
(5) 

{l - 4/i ^F^Vse{t)e{t)^dodoD*{0)} 
(6) 

{l - 4Z2 /, ^F^VMt)e{t)^dodoD*{0)} 

Bringing these contributions in the exponent, we get a standard rotator like 
in eq. Q. Only the inertia has corrections: 

(47) 

+ ^ dodoD*{x)D*{x)D*{x) + -^(8/1 + mh)dodoD*{0) .} 

In the following we consider = 4, corresponding to two-flavor QCD. Since 
D*{0) and dodoD*{0) are known, the only remaining task is to determine the 



'^The NL result is in agreement with that in [J- 
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integral in eq. (|Tf|) . This integral is UV-divergent. The singularity should be 
canceled by the divergent part of the low energy constants h and l2- The 
singularities of these constants are known in dimensional regularization 0]. 



5 Green's functions 



We are in d dimension, d — A -\- e. The physical space-time in Euclidean space 
is Ls X Ls X Ls X [Lt — oo) and the corresponding Green's function is D{x] d). 
Subtracting the ultraviolet and infrared (notation: 'bar' and 'star', respectively) 
divergences we obtain a finite Green's function 

m 

We consider the chiral limit m 0. The UV subtraction (the second term on the 
r.h.s. above) is zero in dimensional regularization . Therefore, L'*(0) = D*{0) 
and 9o9o-D*(0) = 9o9o-D*(0). For the latter, in the framework of DR, the IR 
subtraction is zero also: (?o(^o-C'*(0) — dodoD{0). 

For numerical purposes the following representation is usefujf] 

1 r°° 1 2 f r 21 1 > 

D*ix) = — dw—e--^yo I [n^Siw, y,) - e"™'' \ ^ (49) 

where y = x/Ls and 

oo 

S{w,z) = ^ e-'^»"'cos(27rnz) , 0<w<l, (50) 

n— — OO 

oo 

S{w,z) = l<w<0O. (51) 



n— — oo 



6 The integral J^dx dodoD*{x)D*{x)D%x) 

We calculate this integral using dimensional regularization (DR) in the chiral 
limit. The region of integration is Ls x Lg x Lg x {Lt — >■ oo). We write 

D*{x)^A{x) + D*{x), (52) 

where A(a;) is the infinite volume propagator 

= d,doA{x)^2r^{4xl-r')^^. (53) 

*See, for example, in [14l ITs] . 
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The Green's functions D*{x) and dodoD*{x) are free of IR and UV singularities. 

The integral falls then into six terms. The following four terms (1,2,3,6) are 
finite in DR : 

dx {dodoD*{x)D*{x)D*{x) + dodQA{x)D* {x)D* {x)+ 

^ , (54) 

2A{x)dodoD*{x)D*{x) + dodaA{x)A{x)A{x) \ , 

while the terms below (4,5) are UV singular: 

dx\^A{x)A{x)dodoD*{x) + 2dodoAix)A{x)D*{x)} . (55) 



It is useful to divide the integration region r into a 'cube' and a 'left-right' re- 
gion, where 'cube'= {-Lj2,Ls/2)'^ and 'right'=(-L^/2, L^/2)3 , (L.,Lt) while 
'left'=(-i,/2,i,/2)3 , {~Ls, ~Lt). The left-right region is UV-safc|l. 

Consider the integrals in eq. . The first and third integrals are finite and 
can be integrated as they are. For the second integral we write 



dxdodoA{x)D*{x)D*{x) = 

r r ^ X ^^^^ 

D*{Of J dxdQdoA(x)+ J dxdodoA{x)[D*{xf - D*{Qfj 

The first integral on the r.h.s. of eq. is zero in the cube due to 90° rotation 
symmetry. The rest in eq. (|56|) can be integrated as it is. The fourth integral in 
eq. ([54|) is also zero over the cube. The rest is finite. Let us list the four integral 
values (1,2,3,6) with suppressed: 

-1.228902057d-2, -1.730322906d-2 , 2.461536207d-2 , 4.1588904d-4 . 

(57) 

The final sum of the four integrals in eq. is 

- -^0.00456099852 . (58) 



Turn now to the singular integrals in eq. (|55|) . Write the integral (4) as 

dQdoD*{Q){ [ Aix)A{x)~ [ A{x)A{x) 
^ Jr Jb\s 

+ f A{x)A{x)} + I A{x)A{x){d^d^D*{x)-dodaD*{{))) , 



(59) 



^Tho procedure applied here is similar to that used in |14l I13| . 
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where i? \ 5 is the full space-time with a sphere cut around x = and similarly 
for T \ S. The radius of the sphere is less then, or equal to Lg. Only the first 
integral in eq. (f59|) is singular with DR, the rest is finite. The result is for (4) 



dOdOG* 



1 

8^ 



ln( 



0.009856387107 



(60) 



In the second integrand (5) in eq. ([551) we follow the steps applied above: 
expand D*{x) around x — until the rest gives a finite integral. We write 



D*{x) ^D*{0) + ^\xo\ + {d^d.D*{0)x^x, 



+ {D*ix) D*{0) -1^\xo\-IY1 {d,d^D*iO)x,x,) } 



(61) 



Only the third term on the r.h.s. of eq. (|^T|) is singular in DR. The result is 
1 ( 1 2 



0.015074639535 



(62) 



Collecting the results from eqs. ([55]) . (pO)) . (p^ we obtain 



dxdodQD*{x)D*{x)D*{x) 



(63) 



0.029492025146. 



7 The final result on the inertia O 



The singular part of the bare low-energy constants = 1,2 are known in 
dimensional regularization [4 . For the finite part we follow a generally accepted 
convention |15j . The bare low-energy constants li are written as 



+ I- (64) 



where 7i = 5, 72 = | and 



A=^{^-^(ln(4.)-C-M-hV))}. (65) 
Write the renormalized coupling constants as 
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which gives 



k = 



li 1 7i 



^(ln(47r)-C+l-ln(A2)) 



(67) 



167r2 d~A 327r: 



The combination which enters in eq. (j47p is 



dQdOG* 



2(8^1 + 16^2) = 



dOdOG* r 5 1 



(^^2^2) 



(68) 



-^(ln(4.)-C + l 



) + i(iMA?) + MA^)} 



1 



Multiplying the result in eq. (IM)) by a factor of 4 and adding to eq. (1^ we 
find that the singularities in eq. (|47)) cancel and the result in eq. ([2| is obtained. 
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